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QUANTUM ERGODICITY AND MIXING OF EIGENFUNCTIONS 


STEVE ZELDITCH 


Abstract. This article surveys mathematically rigorous results on quantum ergodic and 
mixing systems, with an emphasis on general results on asymptotics of eigenfunctions of the 
Laplacian on compact Riemannian manifolds. 

Quantum ergodicity and mixing belong to the field of Quantum Chaos, which studies 
quantizations of ‘chaotic’ classical Hamiltonian systems. The basic questions are, how does 
the chaos of the classical dynamics impact on the eigenvalues eigenfunctions of the quantum 
Hamiltonian H and on and long time dynamics generated by HI 

These problems lie at the foundations of the semi-classical limit, i.e. the limit as the 
Planck constant h —0 or the energy E ^ oo. More generally, one could ask what impact 
any dynamical feature of a classical mechanical system (e.g. complete integrability, KAM, 
ergodicity) has on the eigenfunctions and eigenvalues of the quantization. 

Over the last 30 years or so, these questions have been studied rather systematically by 
both mathematicians and physicists. There is an extensive literature comparing classical and 
quantum dynamics of model systems, such as comparing the geodesic flow and wave group 
on a compact (or finite volume) hyperbolic surface, or comparing classical and quantum bil¬ 
liards on the Sinai billiard or the Bunimovich stadium, or comparing the discrete dynamical 
system generated by a hyperbolic torus automorphism and its quantization by the metaplec- 
tic representation. As these models indicate, the basic problems and phenomena are richly 
embodied in simple, low-dimensional examples in much the same way that two-dimensional 
toy statistical mechanical models already illustrate complex problems on phase transitions. 
The principles established for simple models should apply to far more complex systems such 
as atoms and molecules in strong magnetic fields. 

The conjectural picture which has emerged from many computer experiments and heuristic 
arguments on these simple model systems is roughly that there exists a length scale in which 
quantum chaotic systems exhibit universal behavior. At this length scale, the eigenvalues 
resemble eigenvalues of random matrices of large size and the eigenfunctions resemble random 
waves. A small sample of the original physics articles suggesting this picture is P lBU^ lFPl 
IGul UTi IX] . 

This article reviews some of the rigorous mathematical results in quantum chaos, partic¬ 
ularly the rigorous results on eigenfunctions of quantizations of classically ergodic or mixing 
systems. They support the conjectural picture of random waves up to two moments, i.e. on 
the level of means and variances. A few results also exist on higher moments in very special 
cases. But from the mathematical point of view, the conjectural links to random matrices or 
random waves remain very much open at this time. A key difficulty is that the length scale 
on which universal behavior should occur is very far below the resolving power of any known 
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mathematical techniques, even in the simplest model problems. The main evidence for the 
random matrix and random wave connections comes from numerous computer experiments 
of model cases in the physics literature. We will not review numerical results here, but to 
get a well-rounded view of the held it is important to understand the computer experiments 
(see |BSS1 IBarl IKHj for some examples). 

The model quantum systems that have been most intensively studied in mathematical 
quantum chaos are Laplacians or Schrodinger operators on compact (or hnite volume) Rie- 
mannian manifolds, with or without boundary, and quantizations of symplectic maps on 
compact Kahler manifolds. Similar techniques and results apply in both settings, so for the 
sake of coherence we concentrate on the Laplacian on a compact Riemannian manifold with 
‘chaotic’ geodesic how and only briehy allude to the setting of ‘quantum maps’. Additionally, 
two main kinds of methods are in use: (i) methods of semi-classical (or microlocal) analy¬ 
sis, which apply to general Laplacians, and (ii) methods of number theory and automorphic 
forms, which apply to arithmetic models such as arithmetic hyperbolic manifolds or quantum 
cat maps. Arithmetic models are far more ‘explicitly solvable’ than general chaotic systems, 
and the results obtained for them are far sharper than the results of semi-classical analysis. 
This article is primarily devoted to the general results obtained by semi-classical analysis; 
for results in arithmetic quantum chaos, we refer to ini. 

1. Wave group and geodesic elow 

The model quantum Hamiltonians we will discuss are Laplacians A on compact Riemann¬ 
ian manifolds (M, g) (with or without boundary). The classical phase space in this setting is 
the cotangent bundle T*M of M, equipped with its canonical symplectic form ^ 

The metric dehnes the Hamiltonian H{x,^) = where 

~ inverse matrix to [gij\- We denote the volume density of {M,g) 

by dVol and the corresponding inner product on L‘^{M) by (/, g). The unit (co-) ball bundle 
is denoted B*M = {(x,.^) : |.^| < 1}. 

The Hamiltonian flow of H is the geodesic flow. By dehnition, <h*(a:,^) = 
where {xt,^t) is the terminal tangent vector at time t of the unit speed geodesic starting 
at X in the direction Here and below, we often identify T*M with the tangent bundle 
T M using the metric to simplify the geometric description. The geodesic flow preserves the 
energy surfaces {H = E} which are the co-sphere bundles S%M. Due to the homogeneity 
of H, the flow on any energy surface {H = E} is equivalent to that on the co-sphere bundle 
S*M = {H = 1}. (This homogeneity could be broken by adding a potential V G C°^{M) to 
form a semi-classical Schrodinger operator —h^A -|- V, whose underlying Hamiltonian flow 
is generated by |^|g + R(t).) 

The quantization of the Hamiltonian H is the square root \/A of the positive Laplacian, 



of {M,g). Here, g = det[( 7 jj]. We choose to work with \/A rather than A since the former 
generates the wave 

Ut = 
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which is the quantization of the geodesic flow 

By the last statement we mean that Ut is related to in several essentially equivalent 
ways: 


(1) singularities of waves, i.e. solutions t/iV’ of the wave equation, propagate along 
geodesics; 

(2) Ut is a Fourier integral operator (= quantum map) associated to the canonical relation 
dehned by the graph of in T*M x T*M; 

(3) Egorov’s theorem holds. 


We only dehne the latter since it plays an important role in studying eigenfunctions. As 
with any quantum theory, there is an algebra of observables on the Hilbert space L^(M, dvolg) 
which quantizes T*M. Here, dvolg is the volume form of the metric. The algebra is that 
of pseudodifferential operators 'ijjDO^s of all orders, though we often restrict to the 
subalgebra \I/° of 'ipDO^s of order zero. We denote by the subspace of pseudo differen¬ 

tial operators of order m. The algebra is dehned by constructing a quantization Op from an 
algebra of symbols a G S'^{T*M) of order m (polyhomogeneous functions on T*M\0) to T”’’. 
The map Op is not unique. In the reverse direction is the symbol map a a '■ H/”* —^ S'^{T*M) 
which takes an operator Op{a) to the homogeneous term Om of order m in a. For background 
we refer to pCTTllDSjl . 

Egorov’s theorem for the wave group concerns the conjugations 


at {A) ■= UtAU;, A G 


( 1 ) 


Such a conjugation dehnes the quantum evolution of observables in the Heisenberg picture, 
and since the early days of quantum mechanics it was known to correspond to the classical 
evolution 

I/i(a):=ao$* (2) 

of observables a G C^{S*M). Egorov’s theorem is the rigorous version of this correspon¬ 
dence: it says that at dehnes an order-preserving automorphism of i.e. at{A) G 

if A G and that 

■■= Vt{aA), {x,0 e T*M\0. (3) 

This formula is almost universally taken to be the dehnition of quantization of a how or map 
in the physics literature. 

The key difficulty in quantum chaos is that it involves a comparison between long-time 
dynamical properties of and Ut through the symbol map and similar classical limits. The 
classical dynamics dehnes the ‘principal symbol’ behavior of Ut and the ‘error’ UtAU^ — 
Op{aA o ‘h*) typically grows exponentially in time. This is just the hrst example of a ubiq¬ 
uitous ‘exponential barrier’ in the subject. 


2. Eigenvalues and eigenfunctions of A 


The eigenvalue problem on a compact Riemannian manifold 


^jk 

is dual under the Fourier transform to the wave equation. Here, {<Pj} is a choice of orthonor¬ 
mal basis of eigenfunctions, which is not unique if the eigenvalues have multiplicities > 1. 
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The individual eigenfunctions are difficult to study directly, and so one generally forms the 
spectral projections kernel, 

E{\x,y)= ^ ipj{x)ipj{y). (4) 

Semi-classical asymptotics is the study of the \ ^ oo limit of the spectral data {ipj, Aj} or 
of E{X,x,y). The (Schwartz) kernel of the wave group can be represented in terms of the 
spectral data by 

Ut{x,y) = 

3 

or equivalently as the Fourier transform J^e^^^dE{X,x,y) of the spectral projections. Hence 
spectral asymptotics is often studied through the large time behavior of the wave group. 

The link between spectral theory and geometry, and the source of Egorov’s theorem for 
the wave group, is the construction of a parametrix (or WKB formula) for the wave kernel. 
For small times t, the simplest is the Hadamard parametrix, 

poo 

Ut{x,y)^ it<mi{M,9)) (5) 

k=0 

where r{x, y) is the distance between points, Uq{x^ y) = 0“^ (x, y) is the volume 1/2-density, 
inj{M,g) is the injectivity radius, and the higher Hadamard coefficients are obtained by 
solving transport equations along geodesics. The parametrix is asymptotic to the wave 
kernel in the sense of smoothness, i.e. the difference of the two sides of is smooth. The 
relation © may be iterated using Utm = to obtain a parametrix for long times. This is 
obviously complicated and not necessarily the best long time parametrix construction, but 
it illustrates again the difficulty of a long time analysis. 


2.1. Weyl law and local Weyl law. A fundamental and classical result in spectral asymp¬ 
totics is Weyl’s law on counting eigenvalues: 

"W = *{J ■ A, < A} = + 0(A”-‘). (6) 


Here, |H„| is the Euclidean volume of the unit ball and Vol{M,g) is the volume of M with 
respect to the metric g. An equivalent formula which emphasizes the correspondence between 
classical and quantum mechanics is: 


TrEx = 


vomu < A) 

(27r)’^ 


(7) 


where Vol is the symplectic volume measure relative to the natural symplectic form dxjA 

d^j on T*M. Thus, the dimension of the space where H = a/A is < A is asymptotically the 
volume where its symbol |.^|g < A. 

The remainder term in Weyl’s law is sharp on the standard sphere, where all geodesics 
are periodic, but is not sharp on (M, g) for which the set of periodic geodesics has measure 
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zero (Duistermaat-Guillemin inni, Ivrii). When the set of periodic geodesics, has measure 
zero (as is the case for ergodic systems), one has 

A'(A) = #{j ■■ Aj < A} = ^Vol{M,g)X" + o(A”-‘). (8) 

For background, see |TOV] ch. XXIX. The remainder is then of small order than the 
derivative of the principal term, and one then has asymptotics in shorter intervals: 

N([X, A + 1]) = #0 : A, e [A, A + 1]} = n^Vol{M, g)A'‘-‘ + o(A”-‘), (9) 

Physicists tend to write A ~ and to average over intervals of this width. Then mean 
spacing between the eigenvalues in this interval is CnVol{M,g)-^X-^^-^\ where is a 
constant depending on the dimension. 

An important generalization is the local Weyl law concerning the traces TrAE{X) where 
A G It asserts that 

E <Pi) = gW [ cyAdxdOX’' + o(A"-‘). (10) 

A,<A 

There is also a pointwise local Weyl law: 

E (11) 

A,<A ^ 

where R{X,x) = 0{X"'~^) uniformly in x. Again, when the periodic geodesics form a set of 
measure zero in S*M, one could average over the shorter interval [A, A + 1]. Combining the 
Weyl and local Weyl law, we hnd the surface average of cta is a limit of traces: 

j ( 12 ) 

^ ^ X<A 

Here, ja is the Liouville measure on S*M, i.e. the surface measure dfi = induced by 
the Hamiltonian H = |,^|g and by the symplectic volume measure dxd^ on T*M. 

2.2. Problems on asymptotics eigenfunctions. Eigenfunctions arise in quantum me¬ 
chanics as stationary states, i.e. states -0 for which the probability measure |'0(f, a;)p(ino/ is 
constant where 'ip{t,x) = Ut'ip[x) is the evolving state. This follows from the fact that 

Ut^Pk = (13) 

and that = 1. They are the basic modes of the quantum system. One would like to 

know the behavior as Xj —a oo (or h —> 0 in the semi-classical setting) of invariants such as: 

(1) Matrix elements {Aipj,ipj) of observables in this state; 

(2) Transition elements {Aipi, ipj) between states; 

(3) Size properties as measured by norms ||(/3j||Lp; 











STEVE ZELDITCH 


(4) Value distribution as measured by the distribution function Vol{x G M : \(pj{x)\^ > 
t}. 

(5) Shape properties, e.g. distribution of zeros and critical points of (pj. 

Let us introduce some problems which have motivated much of the work in this area. 


Problem 1. Let Q denote the set of ‘quantum limits’, i.e. weak* limit points of the sequence 
{<hfc} of distributions on the classical phase space S*M, defined by 


where a G C°°{S*M). 



{Op{a)ipk,(pk) 


The set Q is independent of the dehnition of Op. It follows almost immediately from 
Egorov’s theorem that Q C Aij, where Aij is the convex set of invariant probability measures 
for the geodesic flow. Furthermore, they are time-reversal invariant, i.e. invariant under 
{x,f) —>• {x, —f) since the eigenfunctions are real-valued. 

To see this, it is helpful to introduce the linear functionals on 

Pk(A) = {Op{a)ipk,H^k)- (14) 

We observe that pk{I) = 1, that pk{A) > 0 if A => 0 and that 

PkiUtAUf) = pk{A). (15) 

Indeed, if A > 0 then A = B*B for some i? G 4/*^ and we can move B* to the right side. 
Similarly is proved by moving Ut to the right side and using (USD. These properties 
mean that pj is an invariant state on the algebra 4/°. More precisely, one should take the 
closure of in the operator norm. An invariant state is the analogue in quantum statistical 
mechanics of an invariant probability measure. 

The next important fact about the states pk is that any weak limit of the sequence {pk} on 

is a probability measure on C{S*M), i.e. a positive linear functional on C{S*M) rather 
than just a state on 4/°. This follows from the fact that {Kipj,ipj) —>■ 0 for any compact 
operator K, and so any limit of {Apki L^k) is equally a limit of ((A -t- K)pki ‘Lk)- Hence any 
limit is bounded by AAk 11 A -|- iL| | (the inhmum taken over compact operators), and for any 
A G 4/°, ||ct^||l°° = infic ||A -|- K\\. Hence any weak limit is bounded by a constant times 
||<ta||l°° and is therefore continuous on C{S*M). It is a positive functional since each pj is 
and hence any limit is a probability measure. By Egorov’s theorem and the invariance of the 
Pk, any limit of Pk{A) is a limit of pk{Op{aA o *h^)) and hence the limit measure is invariant. 

Problem I is thus to identify which invariant measures in A4/ show up as weak limits of the 
functionals pk or equivalently the distributions d<hfc. The weak limits reflect the concentration 
and oscillation properties of eigenfunctions. Here are some possibilities: 

(1) Normalized Liouville measure. In fact, the functional ut of (jl2|l is also a state on 
for the reason explained above. A subsequence {(pj^,} of eigenfunctions is considered 
diffuse if pj,^ —>• eo. 

( 2 ) A periodic orbit measure p.y dehned by p^{A) = -A J^aAds where L.y is the length 
of 7 . A sequence of eigenfunctions for which pkj —/i^ obviously concentrates (or 
strongly ‘scars’) on the closed geodesic. 

(3) A hnite sum of periodic orbit measures. 
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(4) A delta-function along an invariant Lagrangian manifold A C S*M. The associated 
eigenfunctions are viewed as localizing along A. 

(5) A more general measure which is singular with respect to dfx. 


All of these possibilities can and do happen in different examples, 
particular, we have 


If 


uj then in 


Vol{M) 


\(pkj{x)\‘^dVol 


Vol{E) 


Vol{M) 

for any measurable set E whose boundary has measure zero. In the interpretation of 
\ipkj{x)\‘^dVol as the probability density of hnding a particle of energy at x, this says 
that the sequence of probabilities tends to uniform measure. 

However, d^kj —>■ a; is much stronger since it says that the eigenfunctions become diffuse 
on the energy surface S*M and not just on the conhguration space M. As an example, 
consider the flat torus An orthonormal basis of eigenfunctions is furnished by the 

standard exponentials -with k E Obviously, so the eigenfunctions 

are already diffuse in conhguration space. On the other hand, they are far from diffuse 
in phase space, and localize on invariant Lagrange tori in S*M. Indeed, by dehnition of 
pseudodifferential operator, ^e 27 n(fc,a;} _ ^'^^2m(k,x) k) is the complete symbol. 

Thus, 


(Ae- 


27 ri{k,x) 27 ri{k,x)\ _ 

5 ^ / ~ 


a{x, k)dx 


(Xa{x, ^)dx. 
\k\ 


A subsequence of eigenfunctions has a weak limit if and only if tends to 

a limit vector in the unit sphere in M”. In this case, the associated weak* limit is 

(^A{x,^o)dx, i.e. the delta-function on the invariant torus C S*M dehned by 
the constant momentum condition ^ = ^o- The eigenfunctions are said to localize under this 
invariant torus for 

The hat torus is a model of a completely integrable system, on both the classical and 
quantum levels. Another example is that of the standard round sphere S'^. In this case, the 
author and D. Jakobson showed that absolutely any invariant measure u E can arise as a 
weak limit of a sequence of eigenfunctions. This rehects the huge degeneracy (multiplicities) 
of the eigenvalues. 

On the other hand, if the geodesic how is ergodic one would expect the eigenfunctions to 
be dihuse in phase space. In the next section, we will discuss the rigorous results on this 
problem. 

Oh-diagonal matrix elements 


Pjk{A) = (16) 

are also important as transition amplitudes between states. They no longer dehne states 
since Pjk{I) = 0, are no longer positive, and are no longer invariant. Indeed, pjkiUtAU^) = 
so they are eigenvectors of the automorphism Ut of ([Q). A sequence of such 
matrix elements cannot have a weak limit unless the spectral gap \j — Xk tends to a limit 
r G M. In this case, by the same discussion as above, any weak limit of the functionals pjk 
will be an eigenmeasure of the geodesic how which transforms by under the action of 
Examples of such eigenmeasures are orbital Fourier coefficients ^ e~'‘'^^aA(^*(x, ^))dt 
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along a periodic orbit. Here r E We denote by Qr such eigenmeasures of the geodesic 
flow. Problem 1 has the following extension to off-diagonal elements: 


Problem 2. Determine the set Qr of ‘quantum limits’, i.e. weak* limit points of the se¬ 
quence of distributions on the classical phase space S*M, defined by 


where Xj — Xk 



{Op{a)(pk, ^j) 


T -|- o(l) and where a E C‘^{S*M), or equivalently of the functionals pjk- 


As will be discussed in SI the asymptotics of off-diagonal elements depends on the weak 
mixing properties of the geodesic flow and not just its ergodicity. 

Matrix elements of eigenfunctions are quadratic forms. More ‘nonlinear’ problems involve 
the norms or the distribution functions of eigenfunctions. Estimates of the norms 
can be obtained from the local Weyl law m- Since the jump in the the left hand side at A 
is the jump in the right hand side is the jump of R{X,x), this implies 

^ ||^,|U„=0(A“J‘). (17) 

j:\j=X 

For general L^-norms, the following bounds hold on any compact Riemannian manifold |Sog| : 


w^. 


Jllp 


where 


5{p) = 


n 


1 _ ij _ 1 

2 p> 2 ’ 

n-1 (1 _ 1\ 

A2 n'l 


, 2 < p < oo. 

(18) 

2(n+l) ^ 

W <p<^ 

2<p< ^0+1). 

— ^ — n—1 

(19) 


2^2 p' 

These estimates are sharp on the unit sphere S'" C The extremal eigenfunctions 

are the zonal spherical harmonics, which are the L^-normalized spectral projection kernels 
||nX*(^xo)|| centered at any xq. However, they are not sharp for generic (M, g), and it is natural 
to ask how ‘chaotic dynamics’ might influence norms. 


Problem 3. Improve the estimates for {M,g) which ergodic or mixing 

geodesic flow. 


In |SogZ| it is proved that if a sequence of eigenfunctions attains the bounds in (fTTjl . then 
there must exist a point Xq so that a positive measure of geodesics starting at Xq in 
return to xq at a hxed time T. In the real analytic case, all return so xq is a perfect recurrent 
point. In dimension 2, such a perfect recurrent point cannot occur if the geodesic flow is 
ergodic; hence |= o(A^“) on any real analytic surface with ergodic geodesic flow. 
This shows that none of the U’ estimates above the critical index are sharp for real analytic 
surfaces with ergodic geodesic flow, and the problem is the extent to which they can be 
improved. 

The random wave model (see ^ predicts that eigenfunctions of Riemannian manifolds 
with chaotic geodesic flow should have the bounds ||</?a||lp = 0(1) for p < cx) and that 
IIv^aII ]^oo < A/log A. But there are no rigorous estimates at this time close to such predictions. 
The best general estimate to date on negatively curved compact manifolds (which are models 
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of chaotic geodesic flow) is just the logarithmic improvement ||<pjHx,oo = C’(Sx) 
standard remainder term in the local Weyl law. This was known for compact hyperbolic 
manifolds from the Selberg trace formula, and similar estimates hold manifolds without 
conjugate points jBerj. The exponential growth of the geodesic flow again causes a barrier 
in improving the estimate beyond the logarithm. In the analogous setting of quantum 
‘cat maps’, which are models of chaotic classical dynamics, there exist arbitrarily large 
eigenvalues with multiplicities of the order the L°°-norm of the L^-normalized 

projection kernel onto an eigenspace of this multiplicity is of order the square root of the 
multiplicity (jFNDj). This raises doubt that the logarithmic estimate can be improved by 
general dynamical arguments. Further discussion of L‘^-norms, as well as zeros, will be given 
at the end of )|nifor ergodic systems. 

3. Quantum ergodicity 

In this section, we discuss results on the problems stated above when the geodesic flow of 
(M, g) is assumed to be ergodic. Let us recall that this means that Liouville measure is an 
ergodic measure for d)*. This is a spectral property of the operator Id of 0 on L^{S*M, dfi), 
namely that Vj has 1 as an eigenvalue of multiplicity one. That is, the only invariant 
functions (with respect to Liouville measure) are the constant functions. This implies that 
the only invariant sets have Liouville measure 0 or 1 and (Birkhoff’s ergodic theorem) that 
time averages of functions are constant almost everywhere (equal to the space average). 

In this case, there is a general result which originated in the work of A. 1. Schnirelman 

[BEI1IHE2]: 

Theorem 1. jSh.lLIZOl I7t1 IZ.3L irTVl ISnl RTCl IZZwj Let {M,g) be a compact Riemann- 
ian manifold (possibly with boundary), and let be the spectral data of its Laplacian 

A. Then the geodesic flow is ergodic on {S*M,dg) if and only if, for every A G \h°(M), 
we have: 

(1) limA^oo]^EA,<Al(^</^i>V?i) = 0- 

(2) (Ve)(3h)limsup;,^^^^ |(A(pj,(pfc)|2 < e 

This implies that there exists a subsequence of eigenfunctions whose indices jk have 

counting density one for which {Acpj^, ipjfl) —>■ a;(A). We will call the eigenfunctions in such 
a sequence ‘ergodic eigenfunctions’. One can sharpen the results results by averaging over 
eigenvalues in the shorter interval [A, A + 1] rather than in [0, A]. 

There is also an ergodicity result for boundary values of eigenfunctions on domains with 
boundary and with Dirichlet, Neumann or Robin boundary conditions jOLL IHZl IBuj . This 
corresponds to the fact that the billiard map on B*dM is ergodic. 

The first statement (1) is essentially a convexity result. It remains true if one replaces the 
square by any convex function ip on the spectrum of A, 

-uj{A)) (20) 

\j<E 

Before sketching a proof, we point out a somewhat heuristic ‘picture proof’ of the theorem. 
Namely, ergodicity of the geodesic flow is equivalent to the statement that Liouville measure 
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is an extreme point of the compact convex set Aij. In fact, it further implies that oo is an 
extreme point of the compact convex set £’k of invariant states for at of m ; see |Rnj for 
§6.3 for background. But the local Weyl law says that u is also the limit of the convex 
combination Pr extreme point cannot be written as a convex combination 

of other states unless all the states in the combination are equal to it. In our case, uj is only 
a limit of convex combinations so it need not (and does not) equal each term. However, 
almost all terms in the sequence must tend to a;, and that is equivalent to ( 1 ). 

Sketch of Proof of (1) As mentioned above, this is a convexity result and with no additional 
effort we can consider more general sums of the form We then have 

^ (p((A(pfc, ipk) - uj{A)) = ^ - uj{A)ipk, (Pfc)). (21) 

\j<E 

We then apply the Peierls-Bogoliubov inequality 

n 

i=i 

with B = IiE[{A)T — u}{A)]IiE to get: 

V^(((^)t — <v(A) 9 ?fc, ipk)) < Tr <y 9 (n^[(A)r — a;(A)]n^). (22) 

\j<E 

Here, Hg is the spectral projection for H corresponding to the interval [0,h^]. From the 
Berezin inequality [Si, (8.18)] we then have (if <p(0) = 0): 

(1.6.7) ^Ti'v>(nE[(7i)T-w(T]nE) < ^Ti'n£v>([(TT-w(T])nE 

= ^e{,H^{,{A)t — a;(A))). 

As long as ip is smooth, lp[{A)t — <v(A)) is a pseudodifferential operator of order zero with 
principal symbol (p((cr^)T — iv(A)). By the assumption that uje ^ uj we get 

^ ^ E^oo N(E) E p{{ALpk,Pk) - (^{A)) < f p{{aA)T - (^{A))dfi. 

1 ' \.<E J{H=l} 

As T —>■ cx) the right side approaches (p(0) by the dominated convergence theorem and by 
Birkhoff’s ergodic theorem. Since the left hand side is independent of T, this implies that 

^k) - uj{A)) = 0 

for any smooth convex p on Spec(A) with </ 9 ( 0 ) = 0. □ 

As mentioned above, the statement ( 1 ) is equivalent to saying that there is a subsequence 
of counting density one for which > ui. The above proof does not and cannot settle 
the question whether there exist exceptional sparse subsequences of eigenfunctions of density 
zero tending to other invariant measures. To see this, we observe that the proof is so general 
that it applies to seemingly very different situations. In place of the distributions {*hj} 
we may consider the set of periodic orbit measures for a hyperbolic flow on a compact 
manifold X. That is, = ^ /^ / for / G C(W), where 7 is a closed orbit and Ty is 
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its period. According to the Bowen-Margulis equidistribution theorem for closed orbits of 
hyperbolic flows, we have 

where as above /i is the Liouville measure, where is the linear Poincare map and where 
n(T) is the normalizing factor which makes the left side a probability measure, i.e. dehned 
by the integral of 1 against the sum. An exact repetition of the previous argument shows 
that up to a sparse subsequence of y’s, ^ n individually. Yet clearly, the whole sequence 
does not tend to d/i: for instance one could choose the sequence of iterates 7 ^ of a hxed 
closed orbit. 

3.1. Quantum ergodicity in terms of operator time and space averages. The hrst 
part of the result above may be reformulated as a relation between operator time and space 
averages. 

Definition Let A G 6e an observable and define its time average to be: 

(A) := lim ^ f U^AUtdt 

1 —>CXD Z1 J—'J' 

and its space average to be scalar operator 

oj{A) ■ I 

Here, the limit is taken in the weak operator topology (i.e. one matrix element at a time). 
To see what is involved we consider matrix elements with respect to the eigenfunctions. We 
have 

(1.3.2) U;AU,dt^i, v>,) = y (A<p„ <pj) 

from which it is clear that the matrix element tends to zero as T —>■ 00 unless Aj = Xj. 
However, there is no uniformity in the rate at which it goes to zero since the spacing A* — Xj 
could be uncontrollably small. 

In these terms. Theorem □ ( 1 ) says that: 

(A) = uj{A)I + K, where lim ujfiK*K) 0, (23) 

A —>00 

where a;A(A) = TrE{X)A. Thus, the time average equals the space average plus a term K 
which is semi-classically small in the sense that its Hilbert-Schmidt norm square \ \E\K\\\jg 
in the span of the eigenfunctions of eigenvalue < A is o{N{X)). 

This is not exactly equivalent to Theorem^ (1) since it is independent of the choice of 
orthonormal basis, while the previous result depends on the choice of basis. However, when 
all eigenvalues have multiplicity one, then the two are equivalent. To see the equivalence, note 
that (A) commutes with a/A and hence is diagonal in the basis of joint eigenfunctions 
of (A) and of Ut- Hence K is the diagonal matrix with entries {Aipk,g>k) — uj{A). The 
condition is therefore equivalent to 
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Since all the terms are positive, no cancellation is possible and this condition is equivalent 
to the existence of a subset 5 C M of density one such that Qs := {d^k ■ k & S} has only uj 
as a weak* limit point. As above, one says that the sequence of eigenfunctions is ergodic. 

One could take this re-statement of Theorem^ (1) as a semi-classical definition of quantum 
ergodicity. Two natural questions arise. First: 

Problem 4. Suppose the geodesic flow of {M,g) is ergodic on S*M. Is the operator K 
in 

{A)=oj{A) + K 

a compact operator? In this case, a/A is said to be QUE (guantum uniguely ergodic) If 
ergodicity is not sufficient for the QUE property, what extra conditions need to be added? 

Compact would imply that {Kipk,(pk) —^ 0, hence {A(pk,(pk) —*• <v(A) along the entire 
sequence. Quite a lot of attention has been focussed on this problem in the last decade. 
It is probable that ergodicity is not by itself sufficient for the QUE property of general 
quantum ergodic systems. For instance, it is believed that there exist modes of asymptotic 
bouncing ball type which concentrate on the invariant Lagrangian cylinder (with boundary) 
formed by bouncing ball orbits of the Bunimovich stadium (see e.g. |KHj for more on such 
‘scarring’). Further, Faure-Nonnenmacher-de Bievre have shown that QUE does not hold for 
the hyperbolic system defined by a quantum cat map on the torus [FNTTj . Since the methods 
applicable to eigenfunctions of quantum maps and of Laplacians have much in common, this 
negative result shows that there cannot exist a universal structural proof of QUE. 

The principal positive result at this time is the proof by E. Lindenstrauss S of the 
QUE property for the orthonormal basis of Laplace-Hecke eigenfunctions eigenfunctions on 
arithmetic hyperbolic surfaces. It is generally believed that the spectrum of the Laplace 
eigenvalues is of multiplicity one for such surfaces, so this should imply QUE completely for 
these surfaces. Earlier partial results on Hecke eigenfunctions are due to Rudnick-Sarnak 
lESl, Wolpert and others. For more on Hecke eigenfunctions, see |M] . 

So far we have not mentioned Theorem ^ (2). In the next section we will describe a similar 
but more general result for mixing systems and the relevance of (2) will become clear. An 
interesting open problem is the extent to which (2) is actually necessary for the equivalence 
to classical ergodicity. 

Problem 5. Converse QE: What can be said of the classical limit of a guantum ergodic 
system, i.e. a system for which {A) = uj{A) -|- K where K is semiclassically in the sense 
above, or compact?. Is it necessarily ergodic? 

Very little is known on this converse problem at present. It is known that if there exists 
an open set in S*M filled by periodic orbits, then the Laplacian cannot be quantum ergodic 
(see jMOj for recent results and references). But no proof exists at this time that KAM 
systems, which have Cantor-like positive measure invariant sets, are not quantum ergodic. It 
is known that there exist a positive proportion of approximate eigenfunctions (quasi-modes) 
which localize on the invariant tori, but it has not been proved that a positive proportion of 
actual eigenfunctions have this localization property. 
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3.2. Further problems and results on ergodic eigenfunctions. Ergodicity is also 
known to have an impact on the distribution of zeros. The complex zeros in Kahler phase 
spaces of ergodic eigenfunctions of quantum ergodic maps become uniformly distributed 
with respect to the Kahler volume form jNVl ISZj . An interesting problem is whether the 
real analogue is true: 


Problem 6 . Ergodicity and equidistribution of nodal sets. Let C M denote the nodal 
set (zero set) ofipj, and equip it with its hypersurface volume form dH"'~^ induced by g. Let 
{M,g) have ergodic geodesic flow, and suppose that {ipj} is an ergodic sequence of eigenfunc¬ 
tions. Are the following asymptotics valid? 



^Wol{M,g) 



fdVol. 


This is predicted by the random wave model of ^ An equidistribution law for the complex 
zeros is known which gives some evidence for the validity of this limit formula. Let (M, g) be 
a compact real analytic Riemannian manifold and let be the holomorphic extension of the 
real analytic eigenfunction ipj to the complexihcation Me of M (its Grauert tube). Then if 
the geodesic flow is ergodic and if ipj is an ergodic sequence of eigenfunctions, the normalized 
current of integration -^Z^c over the complex zero set of gflj tends weakly to dd\fg\. This 
current is invariant under the geodesic flow and is singular along the zero section. 

Finally, we mention some results on norms of eigenfunctions on arithmetic hyperbolic 
manifolds of dimensions 2 and 3. It is proved in ESI that the joint eigenfunctions of A 
and the Hecke operators on arithmetic hyperbolic surfaces have the upper bound ||v3jj|oo = 
for all j and e > 0 , and the lower bound ||(^jj|oo A c-^log log Xj for some constant 
c > 0 and inflnitely many j. In nisi it is proved that there exists an arithmetic hyperbolic 
manifold and a subsequence ipj,, of eigenfunctions with contradicting the 

random wave model predictions. 


4. Quantum weak mixing 

There are parallel results on quantizations of weak-mixing geodesic flows which are the 
subject of this section. First we recall the classical deflnition: the geodesic flow of {M,g) is 
weak mixing if the operator V) has purely continuous spectrum on the orthogonal complement 
of the constant functions in Lf{S*M,dp). Hence like ergodicity it is a spectral property of 
the geodesic flow. 

We have: 


Theorem 2 . ([Z.3,4]) The geodesic flow 4)* of {M,g) is weak 
conditions (l)-(2) of TheoremU\ hold and additionally, for any A 


(Ve)(35) limsup-^ 

A^oo 




< e 


< A 


mixing if and only if the 
G T°(M), 

(Vr e M) 


The restriction j k is of course redundant unless r = 0, in which case the statement 
coincides with quantum ergodicity. This result follows from the general asymptotic formula. 
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valid for any compact Riemannian manifold {M,g), that 






I sinTr |2 

I Tt 




( 24 ) 


In the case of weak-mixing geodesic flows, the right hand side —>■ 0 as T —>■ cx). As with 
diagonal sums, the sharper result is true where one averages over the short intervals [A, A + 1]. 


4.1. Spectral measures and matrix elements. Theorem |21 is based on expressing the 
spectral measures of the geodesic flow in terms of matrix elements. The main limit formula 
is: 



lim 

A—>c» 


1 


E 




(26) 


where d/io-^ is the spectral measure for the geodesic flow corresponding to the principal 
symbol of A, a a G C°°{S*M,dfi). Recall that the spectral measure of Vt corresponding to 
/ G is the measure d/i/ dehned by 


{Vtf, f)L^{S*M) 


dg,f{T) . 


The limit formula (123) is equivalent to the dual formula (under the Fourier transform) 

^ = ('^tcrA,cTA)L2(s*M)- (26) 

The proof of (0F)|) is to consider, for A G T°, the operator A^A G with At = UtAUt- By 
the local Weyl law, 

A^oo iV I^Aj 

The right side of (123 dehnes a measure dm a on M and (123 says 



dmA{r) 


{VtCA, 0 'a)l2{S*M) 





Since weak mixing systems are ergodic, it is not necessary to average in both indices along 
an ergodic subsequence: 


Xj—^oo 

Dually, one has 


lim (A*A(pj,(pj) = ^e**(^* ^^'>\{Aipi, = {VtaA,(TA)L^s*M)- 

j —irvs f * 


(27) 


f^T+e 


lim 

A,—^oo 


dfi, 


o-A' 


~- 

i : |Ai —Aj—r|<£ ^ ^ 

For QUE systems, these limit formulae are valid for the full sequence of eigenfunctions. 


( 28 ) 
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5. Rate of quantum ergodicity and mixing 

A quantitative refinement of quantum ergodicity is to ask at what rate the sums in The¬ 
orem tend to zero, i.e. to establish a rate of quantum ergodicity. More generally, we 
consider ‘variances’ of matrix elements. For diagonal matrix elements, we dehne: 

^^(A) (29) 

In the off-diagonal case one may view \ as analogous to \ {A(pj, tpj) —uj{A)\^. How¬ 

ever, the sums in (OKI) are double sums while those of (Ei are single. One may also average 
over the shorter intervals [A, A -|- 1]. 


5.1. Quantum chaos conjectures. First, consider off-diagonal matrix elements. One con¬ 
jecture is that it is not necessary to sum in j in (OHjl : each individual term has the asymptotics 
consistent with (j28j) . This is implicitly conjectured by Feingold-Peres in |FPj (11) in the form 




C'^( 


/ Ej Ej ) 


h 


) 


2 np{E) ’ 


(30) 


where C'a(t) = e cr^)dt. In our notation, Xj = h ^Ej and p{E)dE ~ dN{\). 

There are ~ OA”"^ eigenvalues A* in the interval [Aj — r — e, Aj — r -|- e], so (IHn|l says that 
individual terms have the asymptotics of 

On the basis of the analogy between |(A(pi, ^j)\^ and (fj) —ui{A)\'^, it is conjectured 

in |FPj that 


fA(A) ~ 


Ca-uj{a)i{0) 

X'^~^vol{Q) 


The idea is that (p± = ± ^j) have the same matrix element asymptotics as eigenfunc¬ 

tions when Aj — Xj is sufficiently small. But then 2(A(p+, (p_) = {Aipi, (fi) — {Apj, pj) when 
A* = A. Since we are taking a difference, we may replace each matrix element by {Api, pi) 
by {Api, Pi) —uj{A) (and also for pj). The conjecture then assumes that {Api, pi) —uj{A) has 
the same order of magnitude as {Api, pi) — {Apj, pj). Dynamical grounds for this conjecture 
are given in jFFKAMMj . The order of magnitude is predicted by some natural random wave 
models, as discussed below in ^ 


5.2. Rigorous results. At this time, the strongest variance result is an asymptotic for¬ 
mula for the diagonal variance proved by Luo-Sarnak for special Hecke eigenfunctions on 
the quotient H^/S'L(2,Z) of the upper half plane by the modular group jLSl IFaj . Their 
result pertains to holomorphic Hecke eigenforms, but the analogous statement for smooth 
Maass-Hecke eigenfunctions is expected to hold by similar methods, so we state the result as 
a Theorem/Conjecture. Note that H^/S'L(2,Z) is a non-compact hnite area surface whose 
Laplacian A has both a discrete and a continuous spectrum. The discrete Hecke eigenfunc¬ 
tions are joint eigenfunctions of A and the Hecke operators Tp (see jSaj for background). 

Theorem/Conjecture 1. jLSj Let {pk\ denote the orthonormal basis of Hecke eigen¬ 
functions for H^/S'L(2,Z). Then there exists a quadratic form B{f) on C'“(H^/S'L(2, Z)) 
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such that 


N{X) 




,|2 B{fJ) 1 . 
fdVol\ =—— + 0 ( 3 ^. 


When the multiplier / = (^a is itself an eigenfunction, Luo-Sarnak have shown that 

= C'v.^(0)L(^,(Pa) 

where L{X) ^x) is a certain L-function. Thus, the conjectured classical variance is multiplied 
by an arithmetic factor depending on the multiplier. A crucial fact in the proof is that the 
quadratic form B is diagonalized by the (p\. 

The only rigorous result to date which is valid on general Riemannian manifolds with 
hyperbolic geodesic flow is the logarithmic decay jZ6j 

Theorem 3. For any {M,g) with hyperbolic geodesic flow, 

Xj<X 


The logarithm as usual reflects the exponential blow up in time of remainder estimates for 
traces involving the wave group. It is rather doubtful that such a result is sharp. However, 
in the case of two-dimensional quantum cat maps, with eigenspaces of multiplicty A/ log A, 
there may exist orthonormal bases with rather large rates of ergodicity. 


6. Random waves and orthonormal bases 


We have mentioned that the random wave model provides a kind of guideline for what to 
conjecture about eigenfunctions of quantum chaotic system. In this hnal section, we briefly 
discuss random wave models and what they predict. 

By a random wave model one means a probability measure on a space of functions. To deal 
with orthonormal bases rather than individual functions, one puts a probability measure on 
a space of orthonormal bases, i.e. on a unitary group. We denote expected values relative 
to a given probability measure by E. We now consider some specihc Gaussian models and 
what they predict about variances. 

As a model for quantum chaotic eigenfunctions in plane domains, M. V. Berry suggested 
using the Euclidean random wave model at fixed energy jB]. Let £x denote the space of 
(tempered) eigenfunctions of eigenvalue A^ of the Euclidean Laplacian A on M". It is spanned 
by exponentials with k E M"', \ k\ = A. The inhnite dimensional space £x is a unitary 

representation of the Euclidean motion group and carries an invariant inner product. The 
inner product defines an associated Gaussian measure whose covariance kernel Cx{x,y) = 
F,f{x)f{y) is the derivative at A of the spectral function 


E{X,x,y) = (27r) 


-n / ^ E 

^kl<A 


(31) 


Thus, 


Cx{x,y) = —E{X,x,y) = flTr)' 


'ki=A 


eh^-yX)dS = (27r)-"A”-^ / (32) 


'kl=i 
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where dS is the usual surface measure. With this dehnition, C\{x,x) ~ h In order to 
make E(/(a;)^) = 1 consistent with normalized eigenfunctions, we divide by to dehne 

Cx{x,y) = (2%)-^ [ 

One could express the integral as a Bessel function to rewrite this as r(^^) | A|a:— 1 /| | ^ J n -2 (A| 

y\)- 

Wick’s formula in this ensemble gives: 

E(p(x)V(l/)^ = + 2<0A(a;,i/)^]. 

Thus, in dimension n we have: 

^UjV{x)V{y)^{xfLp{yfdxdy-V^] = J^Cx{x,y)^V{x)V{y)dxdy 

~ x--Woi(nr In In ~ vl^f^xdy. 

In the last line, we used the stationary phase asymptotics 

{27r)-^X^-^ [ e'^<"-^’«>d5~C'„(A|x-i/|)-'^cos(|x-i/|A). (33) 

J\i\=i 

Thus, the variances have order in dimension n, consistent with the conjectures in 

[FRIEFKAMM] . 

This model is often used to obtain predictions on eigenfunctions of chaotic systems. By 
construction it is tied to Euclidean geometry and only pertains directly to individual eigen¬ 
functions of a hxed eigenvalue. It is based on the inhnite dimensional multiplicity of eigen¬ 
functions of hxed eigenvalue of the Euclidean Laplacian on There also exist random wave 
models on a curved Riemannian manifold (M, g), which model individual eigenfunctions and 
also random orthonormal bases izaizii. Thus, one can compare the behavior of sums over 
eigenvalues of the orthonormal basis of eigenfunctions of A with that of a random orthonor¬ 
mal basis. Instead of taking Gaussian random combinations of Euclidean plane waves of 
a hxed eigenvalue, one takes Gaussian random combinations X]yA e[A a+i] eigen¬ 

functions of (M, g) with eigenvalues in a short interval in the sense above. Equivalently, one 
takes random combinations with Ylj IgP ~ These random waves are globally adapted 
to {M,g). The statistical results depend on the measure of the set of periodic geodesics of 
(M, g)] thus, as discussed in [KHZj , diherent random wave models make diherent predictions 
about oh-diagonal variances. 

Fix a compact Riemannian manifold {M,g) and partition the spectrum of a/A into the 
intervals Ik = [k, k + V\. Let 11^ = E{k + 1) — E{k) be the kernel of the spectral projections for 
a/A corresponding to the interval Ik- Its kernel 11^(0;, y) is the covariance kernel of Gaussian 
random combinations Yhj Xj&ik is analogous to Cx{x,y) in the Euclidean case; it is 

of course not the derivative dE{\,x,y) but the diherence of the spectral projector over Ik- 
We denote by N{k) the number of eigenvalues in Ik and put Hk = ranflfc (the range of 11^). 
We dehne a random orthonormal basis of Hk by changing the basis of eigenfunctions of 
A in Hk by a random element of the unitary group U{Hk) of the hnite dimensional Hilbert 
space Hk- We then dehne a random orthonormal basis of lF‘{M) by taking the product over 
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all the spectral intervals in our partition. More precisely, we define the infinite dimensional 
unitary group 

U{'^) = n“ ,£/(•«») 

of sequences (f/i, U 2 , ■ ■ ■), with Uk € U(Hk)- We equip U{oo) with the product 

du^ = 

of the unit mass Haar measures duk on U{'Hk)'- We then define a random orthonormal basis 
of L‘^{M) to be obtained by applying a random element U G U{oo) to the orthonormal basis 
$ = {^Pj} of eigenfunctions of 

Assuming the set of periodic geodesics of (M, g) has measure zero, the Weyl remainder 
results (jHl) and strong Szego limit asymptotics of [nniiT7^ give two term asymptotics for 
the traces flfcAflfc, (flfcAflfc)^ for any pseudo differential operator A. Combining the strong 
Szego asymptotics with the arguments of |Z.5j , random orthonormal bases can be proved to 
satisfy the following variance asymptotics: 


CO 


E(E,:a,64 \{AU^vU^,)-u;{A)\^ ~ {u;{A*A) 




(A) 


E(E 


i^j . Aj ,Aj 


sinT{Xi—Xj—T) 

T{Xi-Xj-T) 


\{AUip,,Uip,)\^ 


~ {2 


I sirirT 
I tT 



N{k) 


sinT{Xi — Xj—T) 
T(Xi-Xj-T) 


}{u;{A*A) 
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